Quark- Antiquark and Diquark Condensates in Vacuum in a 2D Two-Flavor 

Gross-Neveu Model* 



Zhou Bang-Rong 

College of Physical Sciences, Graduate School of the Chinese Academy of Sciences, Beijing 100049, China and 
CCAST (World Laboratory), P.O.Box 8730, Beijing 100080, China 

(Dated:) 

The analysis based on the renormalized effective potential indicates that, similar to in the 4D 
two-flavor Nambu-Jona-Lasinio (NJL) model, in a 2D two-flavor Gross-Neveu model, the interplay 
between the quark-antiquark and the diquark condensates in vacuum also depends on Gs/Hs, the 
ratio of the coupling constants in scalar quark-antiquark and scalar diquark channel. Only the pure 
quark-antiquark condensates exist if Gs/Hs > 2/3 which is just the ratio of the color numbers 
of the quarks participating in the diquark and quark-antiquark condensates. The two condensates 
will coexist if < Gs/Hs < 2/3. However, different from the 4D NJL model, the pure diquark 
condensates arise only at Gs/Hs = and are not in a possibly finite region of Gs/Hs below 2/3. 

PACS numbers: 12.38Aw; 12.38.Lg; 12.10.Dm; 11.30.Er 

Keywords: 2D Gross-Neveu model, quark-antiquark and diquark condensates, effective potential, renormal- 
ization, mean field approximation 



I. INTRODUCTION 



It is well known that the four-fcrmion interaction mod- 
els including the Nambu-Jona-Lasinio (NJL) model [lj in 
4D space-time and the Gross-Neveu (GN) models [2| in 
2D and 3D space-time are useful means to explore dy- 
namical symmetry breaking. The NJL model may be 
used in research not only on chiral symmetry breaking in 
Quantum Chromodynamics (QCD) which is connected to 
the quark-antiquark condensates [11, H, EL IS H] but also 
on color superconducting at low temperature and mod- 
erate baryonic density which is related to the diquark 
condensates Q. In fact, owing to that the four-fermion 
interactions of (qq) 2 -form and (qq) 2 -form may be related 
to each other via the Fierz transformation, in any four- 
fermion model, the above two forms of interactions are 
always allowed to exist simultaneously. As a result, one 
will face to interplay between the quark-antiquark con- 
densates and the diquark condensates even at tempera- 
ture T = and the quark chemical potential \x = 0, i.e. 
in vacuum. To clarify such interplay theoretically is cer- 
tainly interesting for deeper understanding of feature of 
the four-fermion interaction models. 

In a preceding paper 0], we have researched a 4D 
two-flavor NJL model and proven that the mutual com- 
petition between the quark-antiquark and the diquark 
condensates does occur in vacuum and, depending on 
the ratio of the coupling constants Gs and Hs in scalar 
quark-antiquark and scalar diquark channel, the diquark 
condensates could either exist alone, or coexist with 
the quark-antiquark condensates, or not be formed com- 
pletely. Then a natural question is that whether similar 
conclusions can be reached from the similar four-fcrmion 
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interaction models in other space-time dimensions? i.e. 
whether these conclusions are of generality for this kind 
of models? To answer this question, we will generalize 
the discussions in Ref.Q to a 2D GN model in this pa- 
per. For convenience of making a comparison with the 
results in the 4D NJL model given in Ref. [|| , we will take 
the fermions (still called quarks) with zero bare-masses, 
two flavors and three colors and assume that the diquark 
condensates will come from the attractive interactions in 
the color anti-triplet channel. 

The paper is arranged as follows. In Sect[TT]we will give 
the 2D two-flavor GN model with diquark interactions, 
analyze its discrete symmetries and possible spontaneous 
breaking induced by the condensates. In Sect JIIII we will 
derive the renormalized effective potential in the space- 
time dimension regularization approach. In Sect HVl by 
means of the effective potential, the ground states of the 
model will be determined and the mutual competition be- 
tween the quark-antiquark and the diquark condensates 
will be expounded. Finally in Sect [V] our conclusions fol- 
low. 



II. 2D TWO-FLAVOR GROSS-NEVEU MODEL 
WITH DIQUARK INTERACTIONS 

The Lagrangian of the model is expressed by 

C = qi^d^q + G s [(qq) 2 + (qij 5 fqf} 

+Hs(qn 5 Ts^Aq c )(q c h5Ts^Aq), (1) 

where, in 2D space-time, q = q a k{t, x) will be two compo- 
nent spinor fields with the two flavors a — u,d and three 
colors k — r (red), g (green), b (blue) and the matrices 
7 p (/i = 0, 1), 75 and C will be 2 x 2 matrices which are 
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defined by 
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and obey the relations 

{/,7 , '}=2r, {75,7^1 = 0, C7^ T C = 7 M . (3) 

The matrix C is related to the charge conjugates of the 
quark fields q by 



q C = Cf, q C = q T C. 



(4) 



The Pauli matrices r = (ti,T2,T3) act in the two-flavor 
space and t$ = (tq = 1,ti,t$) are the flavor-triplet 
symmetric matrices. The \a(A = 2,5,7) are the color- 
triplet antisymmetric Gell-Mann matrices acting in the 
three-color space. It is indicated that the product ma- 
trix C^tsXa is antisymmetric and the term with Hs 
in Eq.(l) corresponds to the diquark interactions which 
could lead to scalar color-antitriplet diquark condensates. 

Since a continuous symmetry can never be sponta- 
neously broken in 2D space-time by the Mermin-Wagner- 
Coleman theorem Q, we are interested only in the dis- 
crete symmetries of the model and their possible sponta- 
neous breaking. In this respect, it is easy to check that 
the Lagrangian (1) has the the following discrete sym- 
metries: R : q(t,x) — q(t,x), parity V : q(t,x) —> 
~f°q(t, —x), time reversal T : q(t, x) -5- j°q(~t, x), charge 
conjugate C : q(t,x) — > q c {t,x), discrete chiral symme- 
try xd ■ q(t, x) ^5 7 5 g(t, x), special parity Vi : q(t, x) 
7 1 g(t, —x), the center Z% of SU C (3): 



III. RENORMALIZATION AND EFFECTIVE 
POTENTIAL 

The coupling constants Gs and Hs in Eq.(l) are di- 
mensionless, thus the theory is perturbatively renormal- 
izable. We will conduct the renormalization operation in 
the dimension regularization approach [ljj. To this end, 
on the basis of the standard procedure, we change the 
dimension of space-time from 2 to 

D = 2 - 2e (7) 

and write the renormalized Lagrangian by 

C D = qi^d t ,q + G S M 2 - D Z G [(qq) 2 + (qi l5 Tq) 2 } 

+H s M 2 - D Z H {q C l l5 T S \ A q){qi 1 r T S \ A q C )X$) 

where M is a scale parameter with mass dimension and 
the 7 M are now 2 D / 2 x 2 D I 2 matrices. The Zq and Z H are 
renormalization constants which will be appropriately se- 
lected to cancel the ultraviolet (UV) divergences in the 
loop integrations. Our main interest lies in interplay be- 
tween the two condensates cf> and <5, so we will derive the 
renormalized effective potential of the model containing 
the order parameters from the two condensates. Define 
the order parameters 

a = -2GsM 2 - D Z G <j), A = -2H S M 2 - D Z H S, (9) 

then in the mean field approximation 11] which is equiv- 
alent to the linearization of the interaction terms in the 
existence of the corresponding condensates [11] i-e. 

(qq) 2 ~ 2<j>qq - <f 
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The subscripts c and / indicate that the corresponding 
matrices act in flavor and color space respectively. 

Assume that the four-fermion interactions could lead 
to the scalar quark-antiquark condensates 



(qq) = <t> 



(•5) 



which will break the discrete symmetries xd and Pi, and 
the scalar color-antitriplet diquark condensates 



(q C ijr l f X 2 q) = S 



(6) 



which will break, besides xd and V\, the symmetry 
down to Z 2 corresponding to r and g color degree of free- 
dom. For the condensates 5, we have made appropriate 
transformations respectively in flavor and color space and 
rotated it into t$ = tq = 1 and A = 2 direction. 



(q ty ij5^2q)(qi"(5^2q c ) — S(qij 5 \2q L ' )+(q c ' ij5^q)S* -\S\ 2 

and in the Nambu-Gorkov basis [131 ] with the bispinor 
fields 



(10) 



we can write 

C D = ^>(x)S- 1 (x)^(x)- 



4G S M 2 - D Z G \H S M 2 - D Z H 



In the momentum space, the inverse propagator S (x) 
for the quark fields will be changed into 



V -i-f 5 X 2 A* i>-a 



7%- (11) 



The effective potential in vacuum corresponding to Cd 
can be expressed by 

_2 I A 12 

F(<t,|A|) 



4G S M 2 - D Z G 4H S M 2 ' D Z H 



+1 



J ^nl^^s- 1 (p)s (p), (12) 
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where 



S (p) = 




represents the propagator for massless quark fields (up to 
a factor i) and the Tr is taken over flavor, color, Nambu- 



Gorkov and Dirac spin degree of freedom. 

It is seen from the structure of the matrix X 2 that the 
blue quarks do not participate in the diquark condensates 
and irrelevant to the order parameter A, thus we may 
conduct similar derivation to the one made in Ref. [8| 
and obtain 



( IJ 4GsM2-dz g + 4H s AP-°Z h + J (2ir)» \ p 2 + is + p 2 + is ) ' 1 j 

I 



where in the third term of the right-handed side, the fac- then it is easy to calculate 
tor 4 come from (the flavor number) x (the color number) 
of the red and green quarks participating in the diquark 

condensates, the factor 2 from the flavor number of the . 2 \ f a ^ 2 dl(u ) f a . 2 f d P 1 

blue quarks and the factor 2 D I 2 is due to the represen- ' J Q du 2 J U J (2tt) d p 2 + u 2 

tation's dimension of the 7^ matrix. For the momentum ^, D/2) (a 2 ) D / 2 

integrations, we may make the Wick rotation, use the = — j—-- — (15) 

Euclidean variable p° = ip°, p l = p l (i = 1, • • • , D - 1) {Atz) d / 2 D/2 
and define the integration 

/jD- -2,2 By means of Eqs. (14), (15) and (7) we can transform 

— ^ln ^ 9 ° with I{a 2 = 0) = 0, (14) Eq.(13) to 
(27r)^ p z 
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where 7 = 0.5772 is the Euler constant. To eliminate 
the UV divergent term with the pole 1 je (up to one- loop 
order), we use the minimal substraction scheme [Toj and 
define the renormalization constants Zq and by 



Z G = 1 



6G S 1 



Z H = 1- 



AH s l 



w e 



(17) 



(18) 



In this way, the renormalized effective potential up to 
one-loop order (after the limit D — > 2 is taken) becomes 



V(a,\A\) 



f-[21n. 



M 2 



AGs 2n 
|A| 2 |A| 2 



In- 



■|A| 2 
M 2 



In 



M 2 



where the denotation 

M 2 = 2ixe^M 2 



•3 
(19) 

(20) 



has been used. 



IV. GROUND STATES 

The effective potential V(a, |A|) depends on the two 
real order parameters a and |A|. Its minimum points 
can be found out analytically, then the ground states of 
the model will be determined. From the extreme value 
conditions dV(a, | A\)/da = and dV(a,\A\)/d\A\ = 
we get the equations that 



1 



2G S 



21n- 



M 2 



\A\ 2 



In- 



M 2 



and 



1 



AH S 



hi- 



M 2 



\A\ 2 



0. 



(21) 



(22) 



It is seen from Eqs. (21) and (22) that the non-zero so- 
lutions of a and |A| will depend on the coupling con- 
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stants Gs and Hs- Since 2D NG model is an asymptot- 
ically free theory and G5 and Hs can become running, 
hence it seems that the non-zero a and |A| will depend 
on the scale parameter M. However, we will prove in 
Appendix that the non-zero a and |A| are practically 
scale- independent. The same result also arose in the 2D 
GN model without the diquark interactions [l4[ and it 
simply reflects the characteristic of the mean field ap- 



proximation. 

For discussion of the feature of the extreme value 
points of Via, |A|) we define the determinant 



K 



A B 
B C 



(23) 



where 



.4 
B 
C 



d 2 V(a,\A\) 1 ( 4Gc 

da 2 2G S \ tt 

d 2 V(a,\A\) d 2 V(a,\A\) _ 

dod\A\ 
d 2 V(a,\A\) 



In - 
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Eqs.(21) and (22) have the following four different solu- 
tions which will be discussed successively. 
1) (<r, |A|) = (0,0). In this case we can calculate and get 
that 



.4 



B = 0, C = -00 and K = oo 2 > 0, 



hence (0,0) is a maximum point of V{a, |A|). 

2) (cr, |A|)=(eri,0). By means of the equation obeyed by 

non-zero a\ 



1 



6G S . M 2 
In — 7T 



0. 



it is easy to obtain that 



A tt' K (2^5 3G S 



.4. 



Obviously, (<Ti,0) is a minimum point of V{a, |A|) if 
Gs/Hs > 2/3 and it will be neither a maximum nor 
a minimum point of V(a, |A|) if Gs/Hs < 2/3. 
3) (cr, |A|)=(0,Ai). Now by using the equation satisfied 
by non-zero Ai 



1 



we get 



.4 



2G.c 



AH S 



AH S . M 2 . 
In —tt- = 0, 

TT Af 



*-iln^ 



, K = -A. 

TT 



If Gs ^ we will have that A = —00, K = —00 thus 
(0,Ai) is not an extreme value point of V(a, |A|). Only 
if Gs — 0, it is just possible that A > and K > 
thus (OjAi) becomes a minimum point of V{a, |A|). It 
is indicated that when Gs = we will have cr = by 



Eq.(9) and V(a, |A|) will reduce to V(0, |A|) with the 
single order parameter | A | coming from the pure diquark 
interactions. 

4) (cr, I A|)=(er 2 ,A2), where the non-zero u-i and A 2 are 
solutions of the equations 



1 



2GS ^21n^ 



■In-r =0, 



1 



4H S 



In- 



M 



(To 



= 0. 



(25) 



(26) 



"2 "I" ^2 

By means of the two equations we may obtain 



A = - I 1 

TT 



2a 2 



A 2 
^2 



> 0. 



> 



Hence (02 , A2) is certainly a minimum point of V(a, | A|) 
as long as the solution exists. For seeking the condition 
in which the solution (02, A2) emerges, we find out from 
Eqs. (25) and (26) that 



<r 2 = M exp 



tt f 1 1 
2 \H~ S ~G~ S 



(27) 



A\ = M 2 exp(- 



4ffs 



1 — exp 



3tt / J_ 



3G.c 



(28) 

Obviously, the solution (02, A2) exists only if Gs/Hs < 
2/3. 

The total results above can be summarized as follows: 
In the 2D GN model given by the Lagrangian (1), the 
locations of the minimum points of the effective potential 
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V(a, | A|) depend on the ratios of Gs and Hs and will be 
at 



(0, Ai) 
(tr,|A|) = < (a 2 , A 2 ) if 
(oi, 0) 



G s /i/s = 
< G s /Hs < 2/3 (29) 
Gs/Hs > 2/3 



These results are very similar to the ones obtained in 
the 4D two-flavor NJL model with both quark-antiquark 
and diquark interactions Q, except the only difference 
that in present model the pure diquark condensate so- 
lution (0, Ai) arises only if Gs = however in the 
4D model it emerges in a limited region of Gs/Hs: 
< Gs/Hs < f(H s ) < 2/3. The critical value 2/3 
of Gs /Hs below which the diquark condensates emerge 
is also due to the same fact that only the two-colors (red 
and green) of quarks participate in the diquark conden- 
sates but all the three colors of quarks get into the quark- 
antiquark condensates. 



V. CONCLUSIONS 

We have proven by effective potential approach that, 
similar to the case of the 4D two-flavor NJL model, in 
the 2D two-flavor GN model with diquark interactions, 
the interplay between the quark-antiquark and the 
diquark condensates depends on the ratio Gs/Hs of 
the coupling constants in scalar quark-antiquark and 
scalar diquark channel. In the ground state of the 
model, only pure quark-antiquark condensates could 
emerge if Gs/Hs > 2/3 and the two condensates could 
coexist if < Gs/Hs < 2/3. However, the pure diquark 
condensates could be formed only if Gs/Hs = 0, which 
is different from the 4D model where the pure diquark 
condensates could arise in a finite region of Gs/Hs- In 
any way, the critical value 2/3 of Gs/Hs below which 
the diquark condensates could emerge is the same for the 
2D GN model and the 4D NJL model and it represents 
the ratio of the color numbers of the quarks participating 
in the diquark and the quark-antiquark condensates. 
This is probably a common characteristic of this kind 
of two- flavor four-fermion models. It is also indicated 
that if the model is considered as some simulation of 2D 
QCD and the four-fermion interactions are assumed to 
come from the heavy color gluon exchange interactions 
-g{q^\ a q) 2 (/j, = 0,1; a = l,---,8) via the Fierz 
transformation [TTJ] , then one will obtain Gs/Hs = 4/3 
for two-flavor and three-color case. Hence based on the 
above results, similar to the 4D NJL model case, we 
will have only the pure quark-antiquark condensates 
surviving and no diquark condensates could appear in 
the ground state of the 2D two-flavor GN model in 
vacuum. 

APPENDIX 

We will prove that in the mean field approximation, in 
the considered 2D GN model the non-zero order parame- 



ters (72 and A2 satisfying extreme value conditions of the 
effective potential V(a, |A|) are scale- independent. First 
rewrite Eqs.(27) and (28) coming from the extreme value 
equations (25) and (26) by 



M exp 
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\H S ' 


' G s )_ 



A 2 2 = M 2 exp( — ir/AHs) — <y 2 , 



(A.l) 



then we will replace the coupling constants Gs and Hs 
by respective running form. For this end, define the bare 
couplings 

G° s = G S M 2 - D Z G , H° s = HsM 2 ~ D Z H . (A.2) 

When viewed as independent variables, G s and Hg do 
not depend on the scale parameter M, i.e. they obey the 
equations 



dO] L= dHl = 
dM dM 



(A.3) 



By means of Eqs. (2.17) and (2.18) which define Zq and 
Zfj and the standard renormalization group method [lol ] , 
we may derive the equations 



M 



M 



dG s 


= M dGs = 


12 C 2 


dM 


dM 


dH s 


- dH s 






= M - = 




dM 


dM 



(A.4) 



noting that M differs from M only by a constant factor. 
Eq. (A4) will give the running coupling constants 
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= £ln 


M 2 


G S (M) 


M 2 


1 


= im 

7T 


M 2 


H S (M) 


M Q 2 



1 



H S (M ) 



(A.5) 



where Mq is an arbitrarily fixed momentum scale. Now 
in Eq. (AI) we make the replacements 

G s Gs (M), H S -> H S (M), 

which will lead to the corresponding substitutions 

0-2 -> cr 2 (M), A 2 -> A 2 (M). 

Then by using Eq.(A5), it is easy to verify that 



a 2 (M 2 ) = a 2 (M 2 ), A 2 (M 2 ) = A 2 (M 2 ), (A.6) 



i.e. (T2 and A2 are scale-independent indeed. 
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